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Abstract 

Using an attractive Hubbard model we examine spatial variations of superconducting order parameter and local 
charge on a two dimensional lattice. For various band filling we show the effect of destruction of the order parameter 
around a non-magnetic impurity. In case of a half-filled system such destruction is accompanied by appearance of 
characteristic charge variations around the impurity with an isotropic distribution. 
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According to Anderson theorem the effect of non- 
magnetic disorder on s-wave superconductors can be 
neglected unless spatial fluctuations of order param- 
eter are present [1,2]. However this theorem invented 
for conventional superconductors cannot be applied to 
other superconductors with a short coherence length 
[3,4] and for those of anisotropic pairing [5]. Another 
interesting situation, where the influence of disorder 
can be important, has been found in a case with an 
interplay between different long range orders like su- 
perconductivity and charge density wave (CDW). In a 
negative U Hubbard model this happens for half-filling 
[6] . Here finite disorder favours superconducting phase 
against CDW [3,7,8,9]. These results show that the 
effect of disorder on superconductivity can be some- 
times beneficial leading to disorder induced supercon- 
ductivity [9]. It must be also noticed that the charge 
ordering mechanism and its interplay with supercon- 
ductivity are of a great interest itself because of exis- 
tence of charge strips in HTc superconductors [10] . In 
the present note we will examine the effect of a single 
impurity on a superconducting order parameter and 



charge fluctuations around non-magnetic impurity in 
a two dimensional lattice. 

We start from negative U Hubbard hamiltonian 
[3,4,11]: 

(i,j)tT ia 

+ U^nii;nii, (1) 
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where ti denotes an impurity potential located at the 
central site i — Q {ei = eo5o,0 and U on-site attraction 
{U < 0) which is the same at each lattice site i. 

Our actual calculations consist of solving, self- 
consistently, the following Bogoliubov-de Gennes 
equation [12]: 
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where Ei = ei + Uni/2 denotes the renormalized site 
energy. The pairing potential and the local charge 
Tii are to be found self-consistently: 
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Fig. 1. Pairing potential (a) and charge rii {i — < l,j >) 
on a two dimensional lattice around the central impurity 
{U ^ -3i, n ^ 0.85, eo ^ * at (/, j) ^ (0,0) ). 

n. = 2^(|<|V(S'0 + Kl'(l - , (3) 

where v enumerates the solutions of Eq. 2 for a given 
band fiUing n. 

To examine the effect of a single impurity on the 
surrounding lattice sites we have solved the above 
equations (2,3) in the real space using the recursion 
Lanczos algorithms for a superconductor [13]. In Fig. 
1 we show the local distributions of the pairing poten- 
tial Ai (Fig. la) and the charge rii (Fig. lb) around 
impurity located in the center of a 2d lattice for a band 
filling n close but slightly smaller than a half-filed 
situation (n = 0.85). Due to the impurity (eo = t) the 
pairing potential Ai is going down rapidly at the cen- 
tral site (Fig. la) . This change is coupled to variation 
of the local charge rii (Fig. lb). Note that both distri- 
butions of Ai and rii have similar form e.g. they show 
similar anisotropy. The situation is quite different for a 
half filled system n = 1 (Fig. 2). Note that in this case 
the paring potential goes down around impurity more 
smoothly and isotropically (Fig. 2a) in comparison 
to the previous case (Fig. la). Interestingly, destruc- 
tion of pairing is associated with strong oscillation of 
charge Ui (Fig. 2b). It is clear that around the central 
impurity the electron charge form a localized wave, 
with a characteristic size of 8-9 lattice spaces. In case 



Fig. 2. Pairing potential Ai (a) and charge ni {i — < l,j >) 
on a two dimensional lattice around the central impurity 
(U = -3t, n = 1.0, eo = t at = (0, 0) ). 

of finite concentration of impurities we can expect su- 
perconductor with islands of a normal phase [11]. Our 
calculations show the appearance of local CDW. That 
in turn indicates that impurities can contribute to a 
phase separation phenomenon [9]. 
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